ApJ, accepted 17/07/2009 

Preprint typeset using I^T^X style cmulatcapj v. 10/09/06 



A HEURISTIC PREDICTION OF THE COSMIC EVOLUTION OF THE CO-LUMINOSITY FUNCTIONS 
D. Obreschkow, I. Heywood, H.-R. Klockner, and S. Rawlings 

Astrophysics, Department of Physics, University of Oxford, Keble Road, Oxford, OX1 3RH, UK 

ApJ, accepted 17/07/2009 

ABSTRACT 

We predict the emission line luminosity functions (LFs) of the first 10 rotational transitions of 
12q16q m g a i ax j es a t redshift z — to z = 10. This prediction relies on a recently presented 
simulation of the molecular cold gas content in ~ 3 ■ 10 7 evolving galaxies based on the Millennium 
Simulation. We combine this simulation with a model for the conversion between molecular mass 
and CO-line intensities, which incorporates the following mechanisms: (i) molecular gas is heated 
by the CMB, starbursts (SBs), and active galactic nuclei (AGNs); (ii) molecular clouds in dense or 
inclined galaxies can overlap; (iii) compact gas can attain a smooth distribution in the densest part of 
disks; (iv) CO-luminosities scale with metallicity changes between galaxies; (v) CO-luminosities are 
always detected against the CMB. We analyze the relative importance of these effects and predict the 
cosmic evolution of the CO-LFs. The most notable conclusion is that the detection of regular galaxies 
(i.e. no AGN, no massive SB) at high z > 7 in CO-emission will be dramatically hindered by the weak 
contrast against the CMB, in contradiction to earlier claims that CMB-heating will ease the detection 
of high-rcdshift CO. The full simulation of extragalactic CO-lines and the predicted CO-LFs at any 
redshift can be accessed onlin43 and they should be useful for the modeling of CO-line surveys with 
future telescopes, such as ALMA, the LMT, or the SKA. 
Subject headings: 



1. INTRODUCTION 

An increasing body of evidence suggests that molec- 
ular hydrogen (H2) widely dominated over atomic hy- 
drogen (H I) in the regular galaxies of th e early uni- 
verse (e.g. lObreschkow fc Rawlingd [2009bh . Empirical 
corner stones towards this conclusion were the measure- 
me nt of strong CO-lin e emission in distant regular galax- 
ies (jDaddi et al.ll2008l ). the detection of HI via Lyman-q 
absor ption against distant quasars (e.g. iProchaska et al.l 
2005), the observational confirmation of a correlation 
between the interstellar g as pressure and H2/HI-ratios 
(|Blitz fc Ro solowskv 20Qa) , and the observational confir- 
mation that g alaxy sizes increase s ignificantly with cos- 
mic time (e.g. iBouwens et aT1l2004D . 

In light of future millimeter/submillimcter telescopes, 
such as the Atacama Large Millimeter/submillimeter Ar- 
ray (ALMA) or the Large Millimeter Telescope (LMT), 
much attention is directed towards the possibility of de- 
tecting the suspected molecular gas (mostly H2) at high 
redshift via the characteristic emission lines of the CO- 
molecule. However, the case for frequent CO-detections 
in regular high-redshift galaxies is by no means secure, 
since neither the cosmic evolution of the H2-mass func- 
tion (MF), nor the evolution of the relationship between 
H2-masses and CO-line luminosities is well constrained 
to-date. An elucidation of this situation seems neverthe- 
less within reach, owing to a long list of specific discover- 
ies over the past two decades (Section [3]), based on which 
computer simulations could already predict the CO-linc 
emission of individual h i gh-redshift galaxies in some de- 
tail ( e.g. ICombes et al.l 1999t iGreve fc Sommer-Larsen 
120081) Moreover. iBlain et alJ (|2000h and ICarilli et al. 
(2002) predicted the number of detectable CO-sources in 
various frequency ranges. They assumed that the CO- 

1 http://www. g-vo.org/Millennium/, prior registration required 



line luminosities evolve with the far-IR (FIR) luminos- 
ity, and they tackled the cosmic evolution of the FIR- 
luminosity function by considering a pure density evolu- 
tion. While this approach is perhaps justified at low red- 
shifts, it probably oversimplifies the physical complexity 
of CO-emission at high redshift (z > 1) as we shall show 
in this paper. 

A missing jigsaw piece in the bigger picture is a physi- 
cal prediction of the cosmic evolution of the galaxy lumi- 
nosity functions (LFs) for different CO-emission lines. In 
this paper, we will attempt such a prediction by concate- 
nating many specific empirical and theoretical findings 
about H 2 and CO. The two main steps towards our pre- 
diction of the CO-LFs are (i) a model for the cosmic evo- 
lution of the H2-MF and (ii) a model for the conversion 
between H2-masses and CO-line luminosities. In this pa- 
per, we shall focus on the latter, whil e adopting the H2- 
mass es of a sample of ~ 3 TO 7 galaxies ([Obreschkow et aLl 
l2009h . simu lated based on the Millennium dark matter 
simulation (Spring e! et al.ll2005h . 

In Section [21 we summarize the galaxy simulation 
producing the H2-masses and various other galaxy- 
properties related to CO-line emission. Our model for 
CO-linc luminosities is developed in Section [31 Section 
[5] presents the prediction of the cosmic evolution of the 
CO-LFs and discusses their dependence on the mecha- 
nisms listed at the beginning of Section 03 Section [6] 
ranks the relative importance of these mechanisms and 
discusses the limitations of their implementation. A brief 
summary is given in Section [7] 

2. SIMULATION OF THE H 2 -MF 

This se ction summarizes t h e cold gas simulation pre- 
sented in Obresch kow et all (|2009l ). Main results and 
limita tions were discussed in detail by lObreschkow et all 
(|2009f l and lObreschkow fc Rawlmgsl (2009b). 
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The simulation has three consecutive layers. The 
first layer is the Millennium Simulation (|Springel et al.1 

2005) , an iV-body dark matter simulation in a peri- 
odic box of comoving volume (500 hr 1 Mpc) 3 , where 
H = 100 h km s" 1 Mpc" 1 and h = 0.73. The sec- 
ond simulation layer uses the evolving mass skeleton of 
the Millennium Simulation to tackle the formation and 
cosmic evolution of galaxies in a semi-analyt ic fashion 
([Croton et all 120061 : iDelucia fe Blai^oll2007r i. This is 
a global approach, where galaxies are represented by a 
list of global properties, such as position, velocity, and 
total masses of gas, stars, and black holes. These prop- 
erties were evolved using simplistic formulae for mecha- 
nisms, such as gas cooling, reionization, star formation, 
gas heating by supernovae, starbursts, black hole accre- 
tion, black hole coalescence, and the formation of stel- 
lar bulges via disk instabilities. The resulting virtual 
galaxy catalog (hereafter the "DeLucia-catalog" ) con- 
tains the positions, velocities, merger histories, and in- 
trinsic properties of ~ 3 ■ 10 7 galaxies at 64 cosmic time 
steps. At redshift z = 0, galaxies as low in mass as 
the Small Magellanic Cloud are resolved. The free pa- 
rameters in the semi-analytic model were tuned to vari- 
ous o bservations in the local universe (see ICroton et al.l 

2006) . Therefore, despite the simplistic implementation 
and the possible incompleteness of this model, the simu- 
lated galaxies nonetheless provide a good fit to the joint 
luminosity/colour / m orphology distribution of observed 
low-redshift ga laxies (|Cole et al.ll2001l : lHuang et al.ll2003l : 
iNorberg etaLl 12002ft the b ulge-to-black hole mass re- 
lation (lHaring fc Rixl [2004). the Tully-Fisher relation 
(|Giovanelli et al.lll997t ). and the cold gas meta llicity as a 
function of stellar mass (jTremonti et al.ll2004h . 

The cold gas masses of the simulated galaxies are 
the net result of gas accretion by cooling from a hot 
halo (dominant mode) and galaxy mergers, gas losses by 
star formation and feedback from supernovae, and cool- 
ing flow suppression by feedback from accreting black 
holes. The DeLucia-catalog does not distinguish between 
molecular and atomic cold gas, but simplistically treats 
all cold gas as a single ph ase. Therefore, the third simu- 
lation layer, explained bv lObreschkow et alj ((2009), con- 
sists of post-processing the DeLucia-catalog to split the 
cold gas masses of each galaxy into HI, H2, and He. 
Our model for this subdivision mainly relies on three 
empirical findings : (i) Most cold gas in regular spira l 
([Lerov et all l2008| ) and elliptical galaxies ([Yound l2002j ) 
in the local universe resides in flat disks, and there is 
evide nce that this featu re extends to higher redshifts 
(e.g. iTacconi et al.l I2006T ). (ii) The surface density of 
the total hydrogen component (HI+H2 ) is approximately 
described by an axially symmetric exponential profile 
(jLerov et al.ll2008ft . 

£ H (r) = S H exp(-r/r disk ) , (1) 
where r d i sk is the exponential scale length and the 
normalization factor £jj c & n be calculated as Eh = 
M H /(27T7j isk ), where M H is the total mass of cold hy- 
drogen in the disk. (hi) The local H2 /HI- mass ra- 
tio closely follows the gas pressure of the interstellar 
medium outside molecular clouds over at least four orders 
of magnitude in pressure and for various g alaxy types 
(|Blitz fc Rosolowskyll2006t iLerov et al.|l2008h . Based on 
those findings, we ( Obresch kow et al.l l2009t ) derived an 



analytic expression of the HI- and H2-surface density pro- 
files, 



S H i(r) = 
S H2 (r) = 



S H cxp(-r/r disk ) 
1 + ^moi cxp(-1.6 r/r disk ) 
s h -Rmoi exp(-2.6r/r d j S k) 
1 + ^moi exp(-1.6 r/r disk ) 



(2) 
(3) 



where -R^oi ^ s H2/HI-mass ratio at the galaxy center. 
This model was applied to the galaxies in the DeLucia- 
catalog to characterize their HI and H2 content (masses, 
diameters, and circular velocities). The resulting hydro- 
gen simulation successfully reproduces many local obser- 
vations of HI and H2, such as MFs, mass-diameter re- 
lation s, and mass-velocity relations (|Obreschkow et all 
2009). This success is quite surprising, since our model 
for HI and H 2 only introduced one additional free param- 
eter to match the observed average space density of cold 
gas in the local universe (fObreschk ow et al.ll2009f h A 
key prediction of this simulation is that the H2/HI-ratio 
of most regular galaxies increases dramatically with red- 
shift, hence causing a clear signature of cosmic "down - 
sizing" in the H 2 -MF (jObreschkow fc Rawlingsl [2009b). 
i.e. a negative shift in the mass scale with cosmic time. 

The simulated H2-MF at z — approximately matches 
the local H 2 -MF inferred from the local CO (1-0) - 
LF (jKeres et al.l 120031 : lObreschkow fc Rawlingsl l2009d i . 
and the few measurements of CO-line emission from 
regular galaxies at z « 1.5 (jDaddi et alj I2008D are 
consistent with the predicted H2-MF at this redshift 
([Obreschkow fc Rawling s 2009b). Furthermore, the pre- 
dicted comoving space density of H2 evolves proportion- 
ally to the observe d space density of star formation rates 
( e . g . I Hopkins! 2 00 7f ) within a factor 2 out to at least z = 3. 
For those reasons, we expect the simulated H2-MF to 
scale reasonably well with redshift. Yet, at z > 5 the 
simulation becomes very uncertain because the geome- 
tries and matter content of regular galaxies are virtually 
unconstrained from an empirical viewpoint. The young 
age and short merger intervals of these galaxies com- 
pared to their dynamical time scales, may have caused 
them to deviate substantially from the simplistic disk- 
gas model. An extended discussion of these and other 
limita tions at low and h i gh re dshift is given in Section 
6.3 of lObreschkow et al.l (|2009h . 

3. MODEL FOR THE CO/H2 CONVERSION 

Most detections of H2 rely on emission lines originat- 
ing from the relaxation of the rotational J-levels of the 
12 C 16 0-molecule (hereafter "CO"). Appendix [Bl provides 
background information on the inference of H2-masses 
from CO-line measurements and highlights the justifica- 
tion and drawbacks of this method. 

To predict the CO-line luminosities associated with the 
molecular gas masses of the simulated galaxies (Section 
[2|) , we shall now introduce a simplistic but physically mo- 
tivated model for the conversion between H2-masses and 
CO-luminosities at any redshift. This model aims to re- 
spect the following theoretical and empirical constraints: 

1. The temperature of molecular gas depends on the 
temperature of the CMB and on the radiative feed- 
back from starbursts (SBs) and active galactic nu- 
clei (AGNs). 
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2. Molecular clumps can shield each other if they over- 
lap along the line-of-sight and in velocity space. 
This effect may not be negligible in the dense galax- 
ies at high redshift, especially if observed edge-on. 

3. While locally observed molecular gas is organized 
in Giant Molecular Clouds (GMCs), the dense gas 
in compact luminous galaxies, such as Ultra Lumi- 
nous Infrared Galaxies (ULIRGs), is predicted to 
follow a smooth distribution. 

4. The CO-line emission of molecular gas is correlated 
with the CO/H2-mass ratio, i.e. to the metallicity 
of the galaxy. 

5. The CMB presents an observing background. The 
absorption of CO-lines against the CMB may sig- 
nificantly reduce the effectively detectable lumi- 
nosities of CO-emission lines. 

These mechanisms will be modeled one by one over the 
Sections EHHiSl 

3.1. Gas temperature and the CO -ladder 

To model the luminosity-ratios of different CO-lincs, 
we analyzed the CO-spectral energy distributions (SEDs) 
of nine galaxies drawn from the literature (see Fig. [T] 
and references therein). This sample includes local 
regular galaxies, local and distant SBs, and distant 
quasi stellar objects (QSOs). Four of these sources 
(SMM J16359+6612, F10214+4724, APM 08279+5255, 
Cloverleaf H1413+135) are known to be strongly magni- 
fied by gravitational lensing. We assume that this has no 
major effect on the flux-ratios between different CO-lines. 
This assumption relies on the fact that the lensed galax- 
ies are far-infrared-bright objects, which makes it likely 
that the strongly lensed regions include the star-bursting 
ones. Those are also the high-excitation regions, which 
seem to dominate the CO-emission of most CO-lines (see 
discussion of M 82 in this Section) . 

Surprisingly, all nine CO-SEDs are well fitted by a 
model for a single gas component in local thermodynamic 
equilibrium (LTE). In Appendix [Cl we show that the 
frequency-integrated line luminosities (= power) of such 
a model scales with the upper level J of the transition as 

J 4 

L 3 oc [1 - exp(-rj)] -— , (4) 

where T e is the excitation temperature, vco = 115 GHz is 
the rest-frame frequency of the CO(l-0)-transition, and 
tj is the optical depth. The latter scales with J as 

r T = 7.2 r c exp (- *E " CO J " ) sinh (h J ), ( 5 ) 
P \ 2k h T c J \ 2k h T c )' V ' 

where r c is a constant. The factor 7.2 in Eq. (O was 
introduced in order for t c to correspond to the optical 
depth of the CO(l-O) line (i.e. T\ — r c ) at the excitation 
temperature T c » 17 K, which is the lowest temperature 
of our model (see end of this Section) . 

If normalized to the CO(l-0)-luminosity L\, the LTE- 
model of Eqs. (J1J [5|) has two free parameters r c and T c . 
In order to apply this model to the velocity-integrated 
fluxes Sj shown in Fi g. [Tj w e use the relation Lj oc J Sj 
(see Appendix El Eq. lAl5j) . 



We first fitted the LTE-model to the observed CO- 
SEDs individually via x 2 - m i n i m i za tion. The resulting 
1-cr confidence intervals of the temperatures T c equal 
10 — 20% of their best-fit values. Despite this uncertainty, 
a clear dependence of T e on the galaxy types (regular, 
SB, QSO) can de detected (see below). By contrast, the 
parameters r c are poorly constrained. Their best-fit val- 
ues range from 0.5 to 5 with no clear trend amongst the 
different galaxy types, and their confidence intervals are 
such that a single parameter t c for all CO-SEDs seems 
to provide a consistent solution. We therefore tested a 
second model, where all nine SEDs share the same pa- 
rameter r c , and foun d that the Bayesian evidence (e.g. 
ISivia fc Skilling 2006) of this 10-parameter model (1 Xr c , 
9 x T e ) against the 18-parameter model (9 x r c , 9 x T ) 
is "strong" with odds of order 10 5 : 1. We therefore as- 
sume a single parameter r c for all galaxies, emphasizing, 
however, that the actual optical depth tj varies consid- 
erably as a function of J and T e by virtue of Eq. ([5]). 
The best fit to all nine SEDs yields r c = 2, consistent 
with the moderate optical dep ths for different CO-lines 
found bv lBarvainis et all (|1997f ) in the Cloverleaf quasar. 
The excitation temperatures T c corresponding to t c = 2 
are listed in Fig.[TJ The individual reduced % 2 's for each 
galaxy range from 0.5 to 1.2, hence demonstrating that 
the LTE-model with a single parameter r c provides an 
excellent fit to all observed CO-line-ratios. 




Upper J-level of the CO-line 

Fig. 1. — CO-S EDs of nine near by and distant galax- 
ies: APM 08279 + 5255 JWeiss et al.ll2007h . C loverleaf H14 1 3+117 
[IBarvainis et al.l |l997h . BR 1202-0725 l lOmont et al.l IT996I; 
Kawabc et al. 1999), the central region of M82 j Weiss et al. 
l2005bl). NGC253 flGiisten et al.l I20M) , SMM J1635 9+6612 
(Weiss et aLll2005al), IRAS F10214+4724 IIAo et al.1 |200H| 1. IC 342 
llsrael fc Baasi 120031) . M51 JWielebinski et alJ 11999; INieten etTaTI 
1999). Symbols and error bars represent the measurements, but 
to avoid confusion, only some error bars are displayed. The lines 
represent fits of our model for thermalized gas with finite optical 
depth, and the corresponding excitation temperatures T e are indi- 
cated for each galaxy. Different colors are used for regular galaxies 
(blue), SBs (red), QSOs (green). 

This conclusion justifies the use of the single compo- 
nent LTE-model as a working model, but it does not 
imply that this model describes the physical reality of 
molecular gas. In fact, it seems that neither the assump- 
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tion of LTE-conditions, nor that of a single gas compo- 
nent are fully satisfied in reality. 

Firstly, the density of molecular gas is often too low 
to collisionally excite the higher rotational levels to ther- 
mal equilibrium. In this case, the LTE-conditions are 
not met. A more accurate description of the excitation 
state is then provi ded by so-called larg e velocity gradi- 
ent (LVG) models (|de Jong et al.|[l975| ). which are more 
complex than the LTE-model. A plausible explanation 
for the surprising success of the LTE-model is that the 
suppression of high-J emission by sub-thermal excita- 
tion can be approximately mimicked by a thermalized 
gas with a slightly underestimated optical depth, or a 
temperature T c slightly below the kinetic temperature of 
the gas. Another explanation is that in real clouds the 
sub-thermal excitation of high-J states could be com- 
pensated by a minor fraction of much warmer (~100K) 
and denser molecular gas, such as is seen next to the star- 
forming cores in nearby mole cular clouds (e.g. the "ridge" 
in the Orion molecular cloud. [Lis fe Sc hilkc 2003|). In any 
case, we stress that the temperatures T c identified by our 
LTE-analysis should not be considered as very accurate. 
Better models, albeit more complex, can be found in the 
references of Fig. [TJ 

Secondly, the assumption of a single component seems 
to work because in most galaxies one component widely 
dominates the total CO-SED. A good example to illus- 
trate this conclusion and its limitations is the nearby 
starburst M82, for which the CO-SED u p to the CO (7 - 
6)-transition has been presented bv I Weiss et al.l (|2005bf ). 
The CO-SED of the center of M 82 is displayed in Fig. Q] 
and is reasonably well described by a LTE-model (re- 
d uced x 2 = 1.1). Ye t, an in-depth LVG-analysis (Fig. 7 
in lWeiss et al.ll2005bh revealed that the center of M 82 ex- 
hibits a low-excitation (LE) and a high-excitation (HE) 
component with kinetic temperatures of ~ 50 K (perhaps 
higher for the LE component), consistent with the single 
temperature of the LTE-model of 49 K (see Fig. [T|) . In 
terms of velocity-integrated fluxes, the CO-SED of the 
LE component peaks around the upper level J = 3 — 4, 
while the HE component peaks around J = 6 — 7. How- 
ever, the flux from the HE component completely dom- 
inates the combined CO-SED, such that the latter still 
peaks around J = 6 — 7. This domination of the HE 
component justifies the use of a single component as a 
working model. Moreover, the domination of the HE 
component, in which the excitation conditions are close 
to LTE, adds another reason for the aforementioned suc- 
cess of the LTE-model. The limitations of the single com- 
ponent model become obvious, when considering th e CO- 
SED of the entire galaxy M 82 ([Weiss et al.ll2005bf ). The 
exceptionally strong gas outflows from the star-bursting 
center add an additional LE component, which domi- 
nates the total CO-SED up to the CO(3-2)-transition or 
perhaps the CO(4-3)-transition. The success of the sin- 
gle component model for the other three SBs in Fig. [1] 
suggests that the strongly CO-luminous outflows of M 82 
are rather anomalous. 

Despite the above limitations of the single compo- 
nent LTE-model, we shall use this model for the rest 
of this paper for three reasons: (i) given current com- 
putational resources, applying an LVG-model to ~ 10 9 
galaxies (i.e. up to ~ 3 • 10 7 galaxies per discrete time 



step) is highly impractical; (ii) as demonstrated above 
(e.g. Fig. [1]), the LTE-model is a reasonable working 
model in the sense that it can approximately fit most ob- 
served CO-SEDs; (iii) the differences between the LTE- 
model and the LVG-model are often much smaller than 
the uncertainties associated with other mechanisms, such 
as cloud overlap in high-z galaxies, metallicity, or gas 
heating by SBs (see Section 16. 2[) . 

Fig. [1] demonstrates that the characteristic excitation 
temperatures T e increase from regular galaxies to SBs, 
and more so to QSOs. This supports the interpretation 
of gas-heating by radiation from SBs and AGNs (see also 
observations of | Weiss et aT| [2007 and theoretical work of 
iMalonev fc Blackll988[) ~ On the other hand, T c must also 
depend on th e temperature of the CMB at the redshi ft 
of the source ([Silk fc Spaanslll997t iCombes et al.1119990 . 

We assume that in regular galaxies the molecular gas is 
heated by a constant specific power (i.e. power per unit 
gas mass) , representing the intra-cloud radiative heating 
by massive stars and supernovae associated with regular 
star formation efficiencies. This specific power implies 
a minimal temperature To for the bulk of the molecular 
gas. In addition, the CMB represents a background tem- 
perature of Tcmb(z) = (1 + z) ■ 2.7 K. If the radiative 
heating of molecular gas happens via absorption by op- 
tically thick dust, then the resulting gas temperature or 
the CO-excitation temperature is r* w T 4 + 7cmb(^) 4 i 
as can be seen from combining the Stefan-Boltzmann 
law with the conservation of energy. Following the same 
argument, we can also include the heating of SBs and 
AGNs via 

Tc = T + Tq MB (z) + Tg B + Ta GN , (6) 

where T$b and Tagn are galaxy-dependent parameters 
characterizing the estimated temperatures of the molec- 
ular gas, if respectively SB-feedback or AGN-feedback 
were t he only sources of radiative heating. ICombes et al.l 
(1999) pointed out that, if the radiative transfer is me- 
diated by optically thin dust with an optical depth pro- 
portional to A~ 2 , the exponents in Eq. ^ should be 
increased from 4 to 6. In reality the exponents in Eq. ^ 
are therefore likely to be somewhat higher than 4. Yet, 
Eq. ([6]) only depends on the precise value of the expo- 
nents in the few cases where the highest temperatures on 
the right-hand-side are comparable, while otherwise the 
highest temperature completely dominates T e . 

In the following, we require that the specific radia- 
tion power (oc Tg B ) acquired by the molecular gas from 
SBs, increases proportionally to the surface density of 
the star formation rate (SFR) Esf for small values of 
Esf, while saturating at an upper limit, characterized 
by the temperature Tg^ K . This saturation level encodes 
all possible self-regulation mechanisms, preventing fur- 
ther heating, such as the suppression of star formation 
by photo-dissociation of molecular gas. To parameter- 
ize the efficiency of SB-heating, we define the charac- 
teristic SFR-density £g F , at which the specific radiation 
power reaches 50% of the saturation level. A minimal 
parametrization of these requirements is given by the 
function 

t sb = T sB xi s sf/ (S sf + £ SF ), (7) 
which reduces to the linear relation Tg B « 
T g ™ ax4 Es F /£ SF for S SF < EgF- To compute 
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Ssf = SFR/(7rrg F ) for the galaxies in our simula- 
tion, we approximate the characteristic length tsf with 
the half-mass radius r^ If of molecular gas and we use 
the SFRs comput ed by the semi-analytic model (see 
ICroton eT al. 2006). In this model stars can form via two 
mechanisms: (i) quiescent continual star formation in 
the disk, which depends on the cold gas surface density; 
(ii) star-bursting activity in the bulge, which is driven 
by galaxy mergers. We shall use the combined SFRs of 
both modes to calculate £sf, since, in principle, both 
modes are likely to cause inter-cloud radiative heating, if 
the corresponding SFR densities are high enough, i.e. of 
order S| F . 

In analogy to SBs, we parameterize the heating from 
AGNs via 



Tagn = 7Tgn 4 Mbh/(M bh + M BH ), 



(8) 



where Tj^f*) is the maximal CO-excitation temperature 
that can be achieved by AGN-heating, Mbh is the black 
hole mass accretion rate, and M BH is the critical ac- 
cretion rate, where the specific heating power is half 
the maximum val ue. In the semi-ana lytic model of the 
DeLucia-catalog (|Croton et alj 120061 ). black holes can 
grow via two mechanisms: (i) a quiescent mode, whereby 
black holes continually accrete material from a static hot 
halo; (ii) a merger mode, where the black holes of merg- 
ing galaxies coalesce, while accreting additional mate- 
rial from the cold gas disks. The free parameters in this 
model, were adjusted such that the predicted relation be- 
tween black hol e mass and bul g e mas s matches the local 
observations bv lHaring k Rbd j2004). Since our model 

for CO-heating only depends on Mbh, we have implic- 
itly assumed that all growing black holes have the same 
heating efficiency, independent of their growth mode and 
physical parameters, such as the black hole mass - a sim- 
plistic assumption, which may well require a more careful 
treatment as large samples of CO-detected AGN become 
available. 

To finalize our model, we need to estimate the five pa- 
rameters T , T s ™ ax , T£g*, E§ F , and Mg H . To fix T , 
we consider regular galaxies (no SB, no AGN) in the 
local universe (Tcmb{z = 0) = 2.7 K), where Eq. ([6|) im- 
plies that T c is nearly identical to To. From simultaneous 
CO (2-1) and CO(l- O) detections in 35 regular galaxies 
in the local universe, iBraine et alJ (|1993l ) concluded that 
the ratio between the brightness temperature luminosi- 
ties is L2/ LJ = 0.89 with a scatter of only 0.06. Accord- 
ing to Eq. (|A13|) . this is equivalent to L 2 /Li = 2 3 • 0.89, 
which, by virtue of Eqs. JUEJ, implies a one-to-one cor- 
respondence between T c To and r c . If we impose r c = 2 
(see above), then To ~ 17 K, which roughly agrees with 
the excitation temperatures of the regular galaxies M51 
and IC 342 for the same depth parameter (see Fig. [I} . 
We therefore fix T = 17 K. 

We further set the critical star-formation density to 
Sg F = 500 M© yr _1 kpc -2 , consistent wit h observations 
of the nuclear SBs of M 82 (|de Griisll200"l[ ) and NGC 253 
(|Beck k BeckwitUfl984h . For those galaxies Eq. then 
implies that T g ^ ax4 = T| B /2, where T SB « T e = 40-50 K 



(see Fig.HJ, hence T? B ax = 50— 60 K. We therefore choose 
Tg B ax = 60 K. Our chosen value for E| F also compares 
well to the st ar-formation dens i ty ~ 10 3 M© yr _1 kpc -2 
predicted by pThompso n et al.l (|2005h for the optically 



thick, dense regions of star forming disks. 

For AGN heating, we choose TfijSfc = 150 K, assum- 
ing that the QSO APM 08279+5255 represents an ob- 
ject close to the maximal possible heating. The crit- 
ical black hole accretion rate A/ BH is assumed to be 
M BH = 10M Q yr _1 , consistent with the higher Ed- 
dington accretion r ates in the samp le of 121 radio- loud 
quasars studied bv iBao et alJ ([2008). Assuming a stan- 
dard radiative accretion efficiency of 10%, this value for 
M BH corresponds to a black hole mass of 5 • 10 8 M©, 
which is on the order of a typical progenitor of the su- 
permassive black holes found in the massive galaxies in 
the local universe. 

3.2. Overlap of molecular clumps 

A reason, why CO-radiation can be used as a linear 
tracer of molecular mass in nearby galaxies despite its 
optical thickness is that most lines-of-sight to the molec- 
ular clumps in nearby galaxies do not cross other clumps, 
and hence CO behaves as if it were optically thin (see Ap- 
pendix [El). However, at high redshift, galaxies are denser 
(e.g. iBouwens et al.ll2004l) and ca rry more molecular gas 
( Obreschk ow k Rawlingsl 1200 9b") . and thus the overlap 
(in space and velocity) of molecular clumps may become 
significant. Such overlap will (i) reduce the directly visi- 
ble surface area per unit molecular mass, and (ii) increase 
t he effective optical depth of th e CO-radiation. 

IBallv et all (|1987t see also iGenzel k Stutzkil Il989h 
identified and analyzed more than 100 clumps in the 
Orion molecular cloud. Based on these data, we assume 
that the diameters and masses of clumps are approxi- 
mately r c i ump = lpc and M c i ump = 500 M©. Given a 
total molecular mass Mn 2 , the number of clumps is ap- 
proximated by 

MH2 (9) 



N 



M, 



clump 



and the "fractional filling factor" of a single clump can 
be approximated as 



1 r 2 

clump 

9H 2 r 2 



(10) 



where qn 2 is the axes-ratio of the inclined galaxy-disk 1 
as seen by the observer and rn 2 is the radius of molec- 
ular gas in this disk, which we take as the galactocen- 
tric radius, where £h 2 (r) given in Eq. ([3]) equals 10% of 
the maximal surface density. The factor 0.1 accounts for 
the fact that clouds are only considered to "overlap" , if 
they have comparable radial velocities, as otherwise they 
become mutually transparent to CO-line radiation. The 
value of 0.1 is a rough estimate based o n the velocity dis- 
per sion of each clump (1 — lOkrns -1 , IBallv et aUll987l 
andlMaddalen a~et al.lll986T ) and the fact that the relative 
velocities between two clumps can vary from 10 km s _1 (if 
in the same cloud-complex) up to the circular velocity of 
the disk of several 100 km s -1 (if in different parts of the 
galaxy). In a more accurate model, the value 0.1 would 
have to be altered with the inclination of the galaxy. 

Assuming that the clumps are randomly distributed in 
space and frequency, we find that the fractional volume 

1 The simulated DeLucia-catalog does not provide galaxy orien- 
tations. We therefore assign inclinations randomly between Odeg 
(face-on) and 90 deg (edge-on) according to a sine-distribution. 
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of the position-velocity space covered by the N clumps, 
counting overlapping regions only once, is 



F = 1 - (1 - 



(11) 



F € [0, 1] is here called the "filling factor" (although 
other definitions of this term exist) and its expression of 
Eq. (fTTjl can be derived iteratively by realizing that the 
filling factor of i clumps, i > 1, is «i = Ki-i + «• (1 — ^i-i) 
with K\ = n. Since the summed volume occupied by all 
clumps in the position- velocity space equals N k, we find 
that any line-of-sight crossing at least one clump, must 
on average cross 



D 



N i 



i - (i - K y 



(12) 



clumps, which also overlap in velocity space. If clumps 
do not overlap (i.e. B = 1), the emergent CO-line lu- 
minosities are proportional to the number of clumps N, 
and hence proportional to the molecular mass Mh 2 (see 
Eq. [9]). However, if the clumps overlap (i.e. B > 1), the 
directly visible surface area of the molecular gas is pro- 
portional to Mb_ 2 I B and the optical depth increases from 
rj to Btj. 

3.3. Clumpy and smooth molecular gas 

Measurements of CO-emission lines in distant ULIRGs 
revealed that the use of CO-luminosity-to-H2-mass con- 
version factors known from local galaxies leads to 
H2-masses on the order of or larger than the dy- 
namical masses infe r red from the circular velocitie s 
(iScoville et alJ fl99ll: iSolomon fe Vanden Bout! l200l . 
This contradiction and high-resolution CO-maps of 
ULIRGs led to the new understanding that the densely 
packed GMCs at the center of massive compact galaxies 
are unstable against the tidal shear and therefore disin- 
tegrate into a smooth blend of gas and s tars. Detailed 
observ ations and geometrical models of iDownes et al.l 
(| 19931) and IDownes fe Solomon! (|1998l ) uncovered that 
the smooth gas is about 5-times more CO-luminous per 
unit molecular mass. Multiple line observ ations of the 
two n earby ULIRGs Arp 220 and NGC 6240 ([Greye et all 
120091 and references therein) seem to confirm this model, 
but they also demonstrate that molecular gas in the 
smooth phase may coexist with clumpy gas (i.e. with 
GMCs) in the outskirts of the nuclear disk. 

To account for the possibility of smooth molecular gas, 
we assume that CO-luminosities per unit molecular mass 
scale proportionally to the efficiency 



£ /clumpy ~l~ 5 • (1 ./clumpy); 



(13) 



where / c i U mpy is the H2-mass fraction in the regular 
clumpy phase (i.e. in GMCs) and (1 — /dumpy) is the 
H2-mass fraction in the regular phase. We define the 
transition between the clumpy and the smooth gas phase 
at the H2-surface density threshold S c = 10 3 M Q pc~ 2 , 
which is between the highest H2-densities observed in 
the local univer se (~ 10 2 M Q pc- 2 , e.g. NGC 6946, 
iLerov et al.ll2008t ) and the most extreme H2-surface den- 
sities of ULIRGs (~ 10 4 Mp) pc" 2 , e.g. Arp 220 and NGC 
6240, iGreve et al.ll2009f ). Assuming a thickness of the 
nuclear disk o f a few 10 pc, consistent wi th the nuclear 
disk model of IDownes fc Solomon! (|l998), the adopted 
value of value of S c corresponds to a volume density of 



> 10 3 cm~ 3 . This value falls in betw een the two volume 
densities found bv lGreve et a l. (2009) for the smooth and 
dense gas phases of Arp 220 and NGC 6240. 

In the smooth phase, the rL^/HI-mass ratio is much 
larger than unity, so that the H2-surface density Sh 2 (r) 
can be safely approximated by the total hydrogen density 
En(r) given in Eq. ((T|). The H2-mass fraction in regions 
less dense than S r , can then be calculated as 



/clumpy 



1 



if E H > S c , 
otherwise. 



(14) 



We shall assume that the overlap factor B is calcu- 
lated in the same way for the smooth component as for 
the clumpy one, which corresponds to approximating the 
self-shielding of the smooth region by the self-shielding of 
a densely packed distribution of clumps with the same to- 
tal volume and mass 2 . Fig. [5] shows the simulated global 
fraction of H2-mass in the smooth phase and the frac- 
tion of CO-power from this phase as a function of cosmic 
time. The predicted monotonic increase of both fractions 
with redshift clearly reflects the strong density and size 
evol ution of cold gas disks predicte d by the simulation 
(e.g. [Obreschk ow fc Raw lings 2009a!). 

At z = 2, the ELi-mass fraction in the smooth phase is 
about 0.3% (corresponding to a fractional CO-power of 
^ 1%), roughly consistent with the fact that the space 
density of ULIRGs is ~ 1% of the space density of nor- 
mal galaxies at this redshift (|Daddi et alJ l2008f ) . The 
remaining 99% of CO-power at z = 2 in the simulation 
stems from clumpy gas, i.e. from GMCs. This result 
seems consistent with recent observational evidence that 
star formation in many active star forming galaxies at 
z ss 2 is distributed on sign i ficantly larger scal es than in 
ULIRGs (|Daddi et al.ll2008tlGenzel et al.ll2008D and that 
star formation propert ies in high -z galaxies are similar 
to those in GMCs (e.g. iGaol 12009) . 

At z = 5, the CO-power from the smooth gas phase is 
predicted to make up ~ 10% of the total CO-power of all 
galaxies. One might argue that such small fractions can 
be neglected. However, in Section [5.31 we will show that 
the contribution of CO-radiation from smooth molecular 
gas at z = 5 will change the space density of the brightest 
objects in the CO-LFs by an order of magnitude. 




2 3 
Redshift z 



8 9 10 



Fig. 2. — Cosmic evolution of the H2-mass fraction in the smooth 
gas phase summed over all galaxies in the simulation (solid line), 
and the corresponding fraction of CO-line power (dashed line). 

2 The filling factor F of the smooth component turns out to be 
very close to 1, and hence B ■ps N k. 
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3.4. Metallicity 

Various recent observations revealed significant varia- 
tions of the CO-luminosity-to-H2-mass conversion factor 
within and amo ngst the Milky Way (MW) and several 
nearby galax i es (lArimoto et al.|[l996t iBoselli et al.ll2002l ; 
lWilso'nl[l995t IIsraelll200Ct iPariione et all 1200 lft . In gen- 
eral, the CO-luminosity per unit molecular mass turned 
out to be roughly proportional to the mass-fraction of 
metals (]Obreschkow fc Rowl ings 2009cj). Such a depen- 
dence may naively be expected for radiation emitted by 
a metallic 3 molecule like CO, but considering the op- 
tical thickness of this radiation, one could also conclude 
that the CO-luminosity per unit molecular mass is nearl y 
independent of the metallicity (jKutner fc LeuneH l985). 
However, lower metallicities imply a lower dust-to-gas ra- 
tio and hence a more efficient destruction of CO by ultra 
violet (UV) radiation, which can restore a positive corre- 
lation between metallic ities and the CO-luminosities (see 
iMalonev fc Bl ack 1988, and references therein). 

Based on these quantitative measurements and qualita- 
tive theoretical considerations, we decided to scale the lu- 
minosity of all CO-lines proportionally to Z = Mz/M gas , 
where Mi is the metal mass in cold gas and M gas is the 
total cold gas mass (including He). 

3.5. Effective luminosity against the CMB 

The CMB power per unit frequency has its maximum 
within the frequency band covered by the CO-lines. For 
example, at z = the CMB peaks between the CO(l- 
0)-line and the CO(2-l)-line, and at z = 5 the CMB 
peaks between the CO(8-7)-line and the CO(9-8)-line. 
Therefore, the absorption of CO-lines against the CMB 
may significantly reduce the detectable luminosities of 
CO-emission lines. 

Within our assumption that clouds are in local thermal 
equilibrium, Kirchhoff's law of thermal radiation glob- 
ally predicts that the absorptivity of the clouds equals 
their emissivity. Explicitly, if a cloud emitting ther- 
mal radiation with a power per unit frequency equal to 
e(v) u(l>, T e ), where e(v) G [0,1] is the emissivity at the 
frequency v and u(y, T e ) is the power per unit frequency 
of a black body, then the absorbed CMB power per unit 
frequency equals e(u) u(v, Tqmb)- The effective CO-linc 
luminosity measured against the CMB is the difference 
between the intrinsic luminosity emitted by the source 
and the luminosity absorbed from the CMB. If the in- 
trinsic luminosity of the source is given by Eq. l[%]). the 
effective luminosity against the CMB can be obtained by 
replacing the second factor (i.e. the black body factor) 
in Eq. @ by 



J 1 



J 4 



exp 



(Jl 



p vcqJ 
7cmb(z) 



-1 



(15) 



In particular, this expression ensures that no radiation 
can be detected from molecular gas in thermal equilib- 
rium with the CMB, since £(J, Tcmb(z) 7 z) — 0. Alterna- 
tively, Eq. [15] could also be expressed i n terms of bright- 
ness temperatures (e.g. Eq. (14.46) of iRohlfs fc Wilson! 
12004ft . 

3 Here, all elements other than hydrogen and helium are referred 
to as "metals". 



4. COMBINED MODEL FOR CO-LINE 
LUMINOSITIES 

A priori, we departed from the idea that CO-line fluxes 
Lj scale with the total mass of molecular hydrogen Mh 2 . 
Combining this assumption with the models of Sections 
13. 1113. 5i we heuristically suggest that the CO-line lumi- 
nosity (power) of the transitions J— > J — 1 is given by 



L.i = Mb 



k Z e 
B 



[1 



-Br, 



t(J,T e ,z), (16) 



where fc is an overall normalization factor. Consistent 
with common practice, we shall define k in such a way 
that Lj is the line-power, obtained by integrating the 
power per unit solid angle emitted along the linc-of-sight 
over all directions, even though the emission need not be 
isotropic. In fact, our model is explicitly non-isotropic, 
since the parameter B depends on the inclination of the 
galaxy via the axes-ratio qn 2 . 

To estimate the value of fc, we note that in local regular 
galaxies all molecular gas is in the clumpy phase (i.e. e = 
1), molecular clumps barely overlap (i.e. B = 1), and 
the excitation temperatures are given by T e w Tq, hence 
tj = r c = 2 and £(l,T e ,0) = 2.45. Therefore, Eq. {16]) 
for J — 1 reduces to 



L x = 2.12 kZM H2 



(17) 



According to Eq. (fT7|) . k is proportional to the stan- 
dard CO/H2 conversion factor a\ ( Appendix IE]) : in 
fact, combining Eq. (JT7J) with Eqs. JHH IA13|) . yields 



fc = 11.9 k h fco A Z 



-1 „-3 „,-3 



Adopting the value 



o-i = 4.6M f 7)(Kkms~ 1 pc 2 )~ 1 typical for the MW 
(|Solomon fc Vanden BoutJ [2005) . and a cold gas metal- 
licity of Z = 0.01 - 0.02, yields fc = 5 - 10 • 10" 8 W kg" 1 . 

Here, we shall fix the parameter fc such that Eq. (I16| 
applied to the H2-masses of our simulated galax- 
ies (Section |} at z = reproduces the observed 
CO(l-0)-luminosity f unction (LF ) of t he local uni- 
verse as measured by iKeres et"aH (]2003l ) (see Fig. [3]). 
A x 2_mm i m i za tion for the luminosity range L\ > 
10 5 - 5 Jy kms -1 Mpc 2 , i.e. the range where our simulated 
H2-MF is complete, yields 



k = 8- 10 _8 Wkg _1 , 



(18) 



which is indeed consistent with the aforementioned value 
predicted from the CO/H 2 conversion of the MW. 

5. RESULTS 

We have applied the model of Eq. (|16p to the galax- 
ies of the hydrogen simulation described in Section (2] 
The predicted CO-LFs for the first 10 rotational transi- 
tions in the redshift range z — — 10 are displayed in 
Fig. [3] For consistency with observer's practice, the lumi- 
nosity scales refer to velocity-integrated luminosities L Y , 
as opposed to the frequency-integrated luminosities L (= 
power) used in Section [3] The conversion between those 
luminosities depends on the wavelength of the emission 
line as explained in Appendix lAl 

The good match (reduced x 2 = 0.7) between the sim- 
ulated CO (1-0)-LF at z = and the local CO(1-0)-LF, 
inferred bv lKeres et al.l (|2003l ) from FIR-selected sample 
of IRAS galaxies, is due to our tuning of the constant fc 
and the reasonably accurate H2-mass distribution of our 
hydrogen simulation (Section [2]). 
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Fig. 3. — Predicted cosmic evolution of the CO-LFs in the redshift range z = — 10. The thick black line represents the CO(l-O) 
transition, while increasingly thin and red lines represent the increasing ly higher or der transitions up to CO(10-9). Filled points and 
error bars represent the observed CO(1-0)-LF of the local universe {Kcrcs ct al. 200 3j), The open circ le wi th error bars corresponds to the 
CO(2— 1) density estimate based on two detections in regular galaxies at z m 1.5 by[Daddi ct al. (200f|) (see [Obreschkow & Rawlings1 l2009bl 
for further explanations). The differential space density 0(L V ) is defined as the number of sources per unit comoving volume and unit 
log 10 (L v ) with a velocity-integrated luminosity L v . 



The simulation shows a clear signature of cosmic down- 
sizing from z — 2 to z — for all CO-transitions. This 
feature reflects the predicted d ownsizing of H2-masses 
(| Obreschkow fc Rawling s 2009b) . For the particular case 
of the CO(2-l)-LF at z = 2, the simulation result is 
roughly consistent with the space density (open circle in 
Fig. [3]) inferred from two recent CO(2-l)-e mission mea- 
surem ents in normal galaxies at z — 1.5 by Dad dTet al.l 
(2008). We note that the offset of this empirical data 
po int from our simulation is l arger in Fig.[3jthan in Fig. 1 
of lObreschkow fc Rawlingsl (|2009b[ K In t he latter, we 
have c ompared the empirical data point of iDaddi et al.l 
(|2008f) to our simulated H 2 -MF at z = 2. To this end we 
converted the measured CO(2-l)-luminosities into H2- 
masses using the standard CO/H2 conversi on factor for 
ULIRGs, a 2 = 1 M Q (K km s^pc" 2 )" 1 (|Daddi et all 
2008; see definition of a 2 in Appendix [B]). However, 
the model for the CO/H2 conversion of this paper yields 
higher values of a 2 for regular h i gh-red shift galaxies, such 
as those found by |Daddlet_aL| ([20081). 

For higher order transitions (J > 5), the predicted 
downsizing even extends out to z w 5, due to the strong 
dependence of these transitions on SBs and AGNs (Sec- 
tion [571]). In Fig. [31 the "dominant" transition, i.e. the 



one with the maximal velocity-integrated luminosity per 
unit cosmic volume, is indicated at each redshift. The up- 
per J-level of this transition increases with redshift due 
to the combined radiative heating by SBs and AGNs. We 
shall now analyze the dependence of the CO-LFs on the 
individual mechanisms modeled in Section [3J 

5.1. Effects of radiative heating by SBs and AGNs 

Fig. g] compares the simulated CO(1-0)-LF and CO(6- 
5)-LF at z — and z — 8 to the corresponding LFs, if 
either AGN-heating or SB-heating is suppressed. At low 
redshift, both SBs and AGNs have nearly no observable 
ef fect on the CO(1-0) -LF, consistent with the conclusion 
of iKeres et al.l (|2003f) that only the highest luminosity- 
bin of the measured CO(1-0)-LF could indicate a devia- 
tion from a Schechter-function distribution, perhaps due 
to SBs in the sample. 

By contrast, the CO(6-5)-LF appears to be signifi- 
cantly boosted by AGNs at z = 0. In fact, this simu- 
lated LF deviates from a Schechter function and exhibits 
two "knees" , respectively corresponding to a "normal" 
galaxy population (left knee) and a more luminous pop- 
ulation heated by AGNs (right knee). Since the lumi- 
nous end of the CO(6-5)-LF is entirely dominated by 
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Fig. 4.— Effects of SB- and AGN-heating on the CO(1-0)-LF 
(black) and CO(6-5)-LF (blue) at redshifts z = and z = 8. The 
solid lines represent the CO-LFs of the full model, such as shown 
in Fig. |3l while the other lines represent the cases where either 
SBs (dashed), AGNs (dotted), or both (dash-dotted) were sup- 
pressed in the si mulation. The red d ots with error bars represent 
the local HX-LF (Os hima et alj|20011 ). mapped onto the CO(6-5)- 
luminosity scale as explained in Section |5. II 

AGN-heating, we expect the local space density of the 
most CO(6-5)-luminous objects to match the space den- 
sity of local AGNs. To test the simulation, we there- 
fore overlayed the simulated local CO(6-5)-LF with the 
most recent empirical determination of t he local hard 
2 - 8 keV) X-ray-LF (HX-LF) obtained bv lYencho etai] 
20091) (data points in Fig. gj top). This HX-LF re- 
lies on a galaxy sample studied by the X-ray Chan- 
dra observatory. In order to map the HX-luminosity 
scale onto the CO(6-5)-luminosity scale, we crudely as- 
sumed a proportional relation between the two, tuned 
to the empirical data from the Cloverleaf quasar. We 
evaluated the lensed HX-luminosity (at 2 — 8 keV rest- 
frame) of the Cloverleaf quasar directly from the X-ray 
SEP measured and c orrected for Galactic absorption by 
lOshima et al.l (|200lD . For the cosmology of this paper, 
this HX-luminosity is L H x = (1 ± 0.5) • 10 45 ergs- x . 
On the other hand, the lensed CO(6-5)-linc luminosity 
of the Cloverleaf quasar, interpolated f rom the CO(5- 
4) and CO(7-6) line fluxes presented bv lBarvainis et all 
( 1997) and corrected for the standard cosmology of this 
paper, amounts to L 3=6 w 5 • 10 10 Jy kms -1 Mpc 2 . As- 
suming that the LHx/-kj=6-ratio of the Cloverleaf quasar 
is not affected by differential magnification and that 
it mimics the Lhx /£.i= 6-ratio of local AGNs, the HX- 
LF (jYencho et al.l l2009) transforms into the data points 
shown in Fig.[4j The vertical error bars represent the sta- 
tistical density uncertainties given for the HX-LF, while 
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Fig. 5. — Effects of overlapping clumps and smoothly distributed 
gas on the CO(1-0)-LF at z = 5. The solid line represents the 
GO-LF of the full model, such as shown in Fig. [3] The other 
lines correspond to the suppression of overlap effects (dashed) and 
smoothly distributed gas (dotted). Note the different scale of the 
axes compared to the other figures. 

the horizontal error bars represent the 50% uncertainty 
of -Lrx- The good fit between the space densities of local 
AGNs and those predicted for the luminous CO(6-5)- 
sources supports our prediction. 

At very high redshift (z > 7), where the predicted 
space density of AGNs in the DeL ucia-catalog is ex- 
tremely low (see iCroton et al.l |2006[ ) , SBs become the 
dominant source of CO-heating as shown in Fig. [4] (bot- 
tom). This analysis predicts that SB-heated molecular 
gas disks are the most likely objects to be detected in 
CO-line emission at z > 7. The optimal transitions 
are CO (8-7) and CO(6-5) in terms of velocity-integrated 
luminosities or surface-brightness temperatures, respec- 
tively. 

The effects of gas heating by the CMB will are dis- 
cussed together with the effects of the CMB as an ob- 
serving background in Section 15.51 

5.2. Effects of overlapping molecular gas 

The effect of overlapping clumps (Section l3.2p exhibits 
a modest dependence on the upper J-level of the CO- 
transitions, although minor differences may occur due to 
the dependence of the optical depths on J (see Eq. [5]) . 
Our model predicts that the effect of overlapping clouds 
becomes increasingly important with redshift, as a direct 
consequence of the predicted increase in the surface den- 
sities of galaxies with redshift. Between z = and z — 1, 
the effect is negligible (i.e. < 0.1 dex luminosity change), 
while at z — 5 CO-luminosities are predicted to be re- 
duced by a factor 2-3 due to cloud overlap. The dashed 
line in Fig. [5] illustrates the effect of ignoring the overlap 
of clumps (by forcing £? = l)atz = 5. 

5.3. Effects of smooth molecular gas 

In our model (see Section l3~3")) . we assumed that molec- 
ular gas in very dense galaxy-parts is smoothly dis- 
tributed, rather than organized in gravitationally bound 
GMCs. Within our simplistic treatment (Eqs. [EH fTTjj) . 
this effect is independent of the J-level of the CO- 
transition. The effect of smooth gas becomes increas- 
ingly important with redshift, as a direct consequence of 
the predicted increase in the surface densities of galaxies 
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Fig. 6. — Effects of the cosmic evolution of cold gas metallicity on 
the LFs of CO(l-O) (black) and CO(10-9) (blue) at z = 5. Solid 
lines represent the CO-LFs of the full model, such as shown in 
Fig. \3\ while dashed lines represent the CO-LFs, where the cosmic 
evolution of metals has been suppressed in the simulation. 



with redshift. Between z — and z = 1 the effect is neg- 
ligible (i.e. < 0.1 dex luminosity change), but at z = 5 
its importance is comparable to that of heating by strong 
SBs and massive AGNs. The dotted line in Fig. [5] shows 
the effect of ignoring the possibility of smooth gas (by 
forcing / c ium P y = 1) at z = 5. 

5.4. Effects of metallicity 

Fig. [6] shows the effect of neglecting the cosmic evolu- 
tion of cold metals in galaxies, by illustrating the effects 
at z = 5 of suppressing this evolution. In general, the ef- 
fect of metallicity appears to be relatively weak, since the 
cosmic evolution of the cold gas metallicity from z = 5 
to z = is rel atively weak as dis c ussed in more detail in 
Section 6.3 of lObreschkow et all (|2009t ). 

We also note that the cosmic evolution of the cold gas 
metallicity has a stronger effect on weak CO-sources than 
on the luminous ones. In fact, galaxies more luminous 
than the "knee" of the CO-LFs are nearly unaffected 
by the cosmic evolution of metals at z = 5 compared 
to z = 0. The reason for this feature is that the more 
CO-luminous galaxies are, on average, more massive and 
older, and hence they have already formed the bulk of 
their metals at z > 5. 

5.5. Effects of the CMB 

The CMB as an observing background already becomes 
noticeable at z ~ 1, where its effective reduction of the 
CO(l-0)-luminosities amounts to about 0.1 dex (less for 
higher order transition) according to Eq. (|15| . This effect 
increases steeply with redshift and reaches 1 dex at z » 
5 for the CO (1-0) line, such as shown in Fig. [7l The 
increase of the CO-excitation temperature by the heating 
effect of the CMB appears to be a minor effect, which 
only becomes noticeable around z 5. This effect acts 
against the loss of detectable luminosity by the CMB as 
an observing background by about 10% at z w 5 for the 
CO (1-0) line (slightly more for higher order transitions). 

Our model generally predicts that the effect of the 
CMB as an observing background always dominates the 
opposite effect of the CMB as a source of heating. Hence, 
the combined effect of the CMB always reduces the de- 
tectable luminosities of CO-sources, at all redshifts and 




Fig. 7.— Effects of the CMB on the CO(1-0)-LF at z = 1, z = 5, 
and z = 10. The solid lines represent the CO-LFs corresponding to 
the full model, such as shown in Fig. [3] The other lines represent 
the cases, where either the gas- heating by the CMB (dashed), the 
CMB as an observing background (dotted), or both effects (dash- 
dotted) have been suppressed in the simulation. 

for all rotational transitions. This effect is most pro- 
nounced for lower order transitions, since emission from 
higher order transitions originates mostly from SBs and 
AGNs, whose heating effect can mask the comparatively 
low temperature of the CMB (e.g. Tcmb <C Tsb and 
Tqmb "C ?agn)- For regular galaxies (no SBs, no 
AGNs) , the combined effect of the CMB increases with 
redshift in such a way that these galaxies become vir- 
tually invisible in CO-line emission against the CMB at 
z>7. 

Thi s res ult contradicts the claims of Sil k fc Spaansl 
l|1997f) and iGnedin et all {2001) that the higher excita- 
tion temperatures caused by the warm CMB of the early 
universe will ease the detection of CO-emission lines. 
This conclusion only accounts for gas heating by the 
CMB, but ignores the CMB as an inevitable o bserv ing 
background as already noted bv lCombes et al.l (|1999f ) . 
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6. DISCUSSION 

6.1. Ranking of various mechanisms 

The analysis of Section [5] can be summarized in a 
ranking list of the different mechanisms affecting the 
CO/H2 conversion. This ranking naturally depends on 
the redshift, the considered part of the CO-LF, and the 
J-level of the CO-transition. Here we consider the source 
population around the "knee" of the CO-LFs for the lev- 
els J = 3 — 6 at redshift z ~ 3. This ca se corresponds 
to using the fully funded (|Maiolind 120081 ) ALMA-bands 
3 and 4 for the first ALMA-science goal, i.e. the detec- 
tion of a MW-typ e galaxy in CO-line emission at z » 3 
(|De Breuckl 120051) . For this particular setting the pre- 
dicted ranking (from most important to least important) 
is 

1. Gas heating by AGNs (+) 

2. Gas heating by SBs (+) 

3. Overlap of clumps (— ) 

4. Smooth gas (+) 

5. CMB as an observing background (— ) 

6. Cosmic evolution of the cold gas metallicity (— ) 

7. Gas heating by the CMB (+) 

The signs in parentheses indicate whether the effect in- 
creases (+) or decreases (— ) the CO-line luminosities. 
Although this ranking may change considerably with red- 
shift and with J (e.g. Fig. HJ), the above ranking can be 
considered as a rule of thumb for estimating the relative 
importance of various effects. For example, if a simula- 
tion of CO-LFs includes a model for the smooth gas in 
high-redshift galaxies, then it should also account for the 
heating by SBs and AGNs and self-shielding by overlap- 
ping clumps. 

6.2. Model limitations 

The predictions presented is this paper are approx- 
imate ramifications of a semi-empirical model, which 
potentially suffers from simplifications and uncertain- 
ties on each of the four successive simulation-layers: (i) 
the Millennium dark matter simulation, (ii) the semi- 
analytic galaxy simulation of the DeLucia-catalog, (iii) 
our post-processing to assign extended HI- and H 2 - 
properties to each galaxy, (iv) the model for CO-line 
emission introduced in this paper. It is beyond the scope 
of this paper to discuss the limitations related to the 
simulation-l ayers (i)-(iii) , but ext e nsive discussions wer e 
pro vided by iSpringel et all (120051) , iCroton et all (|2006l ), 
and lObreschkow et alH 2009D . respectively. 

All four simulation-layers were widely constrained 
by a broad variety of observations: (i) the cosmo- 
logical parameters for the Millennium si mulation were 
adopted from 2dFGRS dColless et al J 1 2001 and WMAP 
(jSpergel et all [20031: Bennett et aTf l2003): (ii) the semi- 
analytic recipes are motiva ted by various references 
given in ICroton et al.l (|2006l ) and the free parameters 
were tuned to fit the luminosity /c olour / morphology 
distribution of low-redshift galaxies (|Cole et al.l 120011 : 



iHuang et all 120031: iNorberg et ail 120021). th e bulge-to- 
black hole mass r elation (lHaring fe RixH2004f ). the Tully- 
Fisher relation (|Giovanelli et alj 119971 ) , the cold gas 
meta llicity as a function of stellar mass (|Tremonti et al.l 
l2004f ): (iii) our model to assign HI- and ^-properties 
is motivated by va rious studies of HI and H2 in 
local galaxies (e.g . lYoungj 120021: iBlitz fc Rosolowskvl 
I2006t iLerov et all 120081 : lElmegreenl 1 19931) and the 
free parameter was tuned to th e local space den- 
sity of cold gas (iKeres et al.1 [20031: IZwaan et all 120051 : 
lObreschkow fc Rawhngsll2009cD : (iv) our CO-model was 
constrained as decribed in Section [3] of this paper. Since 
this empirical basis is widely dominated by observations 
in the local universe, we expect our low-redshift predic- 
tions for CO, such as the CO-LFs for higher order tran- 
sitions, to be more accurate than the high-redshift pre- 
dictions. 

With regard to our model for CO-line emission, the 
most reliably modeled effects are those of the CMB, since 
they could be assessed from global considerations, with 
no strong dependence on free parameters (see Sections 
13.51 15. 1|) . Also the effects of cold gas metallicity are 
relatively unproblematic: firstly, these effects are rela- 
tively small (e.g. Fig. [S]); secondly, the correlation be- 
tween metallicity and the CO/H2 conversion is empir- 
ically supported (see Section 13. 4[) : thirdly, the metal- 
licities predicted by the DeLucia-catalog seem reliable 
as they reproduce the mass-metallicity relation inferred 
from 53,000 st ar forming galaxies in the Sloan Digi- 
tal Sky Survey (iDe Luc ia et alj|200i ICroton et~aI1l2006t 



Tremonti et al.ll2004D . 

The effects of inter-cloud heating by high density star 
formation (or SBs) and heating by AGNs are much less 
certain. Our temperature model relies on the CO-SEDs 
of only 7 galaxies (4 SBs and 3 QSOs) with poorly con- 
strained star formation densities and black hole accretion 
rates. It is further possible that the molecular material 
in some of these galaxies is simultaneously heated by 
both a SB and an AGN. The relations of Eqs. flBHEJ) are 
simplistic parameterizations of our limited knowledge on 
gas heating by SBs and AGNs, but they may require 
a revision as larger galaxy samples with simultaneous 
CO-SEDs, SFRs, and black hole accretions rates come 
on line. Given the present-day uncertainties of SB- and 
AGN-heating, the use of the LTE-model for CO-SEDs 
(Section 15. ip seems sufficient. In fact, the offset of the 
LTE- model model from the more complex LVG-models is 
small (|Combes et al.l II999I ) compared to the systematic 
uncertainties of radiative heating by SBs and AGNs. 

Surprisingly, we found that self-shielding by overlap- 
ping clumps is perhaps the most subtle effect to model 
at z > 3, because it seems to be a very significant effect 
(see Fig. [5]) and yet its physical complexity is consid- 
erable. Especially in the case of galaxies with heavily 
overlapping clumps (i.e. B 3> 1), the value of the overlap 
parameter B sensibly depends on the radius and mass of 
molecular clumps. If we also consider that clumps are not 
randomly distributed, but organized in cloud-complexes, 
and that their geometries are far from spherical, the pre- 
dicted CO-line luminosities of galaxies with heavily over- 
lapping clumps could differ from our current prediction 
by nearly an order of magnitude. Similar uncertainties 
should be assumed for the effects of smooth gas in high- 
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redshift galaxies. In fact, the critical surface density S c , 
at which gas transforms from clouds to smooth disks, is 
very uncertain and may vary as a function of the mass 
and size of the galaxy. 

7. CONCLUSION 

We have predicted the cosmic evolution of the galaxy 
LFs for the first 10 rotational transitions of the CO- 
molecule. This prediction relies on a combination of a 
recently presented simulation of H 2 -masses in ~ 3 • 10 7 
evolving galaxies with a model for the conversion be- 
tween H2-masses and CO-line luminosities. The latter 
model accounts for radiative heating by AGNs, SBs, and 
the CMB, for smooth and overlapping gas, for the cosmic 
evolution of metallicity, and for the CMB as an observing 
background. 

The main outcome of this study is two-fold. Firstly, the 
predicted CO-LFs are probably the most robust basis to- 
date towards predicting the CO-line detections of high- 
red shift surveys with future telescopes, such as ALMA 
(see lBlain et"aLll200Q[) . the LMT, or phase-3 of the Square 
Kilometre Array (SKA). Secondly, this study revealed 
that the most serious uncertainties of the CO-LFs at 
high redshifts originate from the poorly understood self- 
shielding of overlapping clouds, from the smooth gas in 
luminous galaxies, and from the heating by SBs and 
AGNs. Hence, any serious progress in predicting the 
CO-LFs must address these mechanisms in more detail. 
By contrast, the widely cited effects of the CMB and the 
cosmic evolution of metallicity seem to be relatively well 
modeled. 

This study makes some explicit predictions, which 
could be tested in future CO-surveys; e.g.: 

1. The CO-LFs associated with the first 10 rota- 
tional transitions should show a strong signature 
of "downsizing" in the redshift range z = — 2. 
Explicitly, the total power of each CO-line per co- 
moving volume increases from z — to z — 2 by 
a factor 2 to more than 10, depending on the CO- 



transition (see Fig. [3]). 

2. On average, the relative CO-line power in higher 
order transitions, i.e. the excitation temperature 
T c , increases monotonically with redshift z. This 
is a consequence of more heating at high z, mainly 
due to SBs and AGNs (see Fig. [3]). 

3. Some CO-LFs (e.g. CO(6-5) and CO(7-6) at z = 0, 
and CO(5-4), CO(6-5) and CO(7-6) at z = 2) sig- 
nificantly deviate from a Schechter function. They 
are predicted to have two "knees" , respectively cor- 
responding to a "normal" galaxy population and 
a more CO-luminous population, where the gas is 
heated mostly by AGNs (e.g. Fig. [4]top). 

4. Out to the most distant galaxies, most of the 
cosmic CO-luminosity is predicted to stem from 
regular clumpy gas, i.e. from GMCs, rather than 
a hypothetical dense phase, which is believed 
to dominate some ULIR Gs dDownes et al.1 HqqI 
iDownes fc Solomon! [l998h . 

5. The CMB will significantly suppress the apparent 
CO-line flux of galaxies at high z (see Fig. [7]). In 
particular, galaxies at z > 7 with no strong source 
of internal heating, such as a SB or an AGN, will 
not be detectable in CO-line emission. 
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APPENDIX 

A. LUMINOSITIES AND FLUXES OF LINES 

This section overviews the concepts used in relation with line fluxes and line luminosities with an emphasis on 
connecting the terminology and units of observers to those of theoreticians. 

Terminology and definitions 

Any continuous isotropic electromagnetic radiation from a point-source is completely characterized by the luminosity 
density (or monochromatic luminosity) l{v), an intrinsic quantity measured in units proportional to lWHz -1 = 1 J. 
The corresponding observable quantity is the flux density (or monochromatic flux) s(v), measured in units proportional 
to 1 W Hz -1 m~ 2 = 1 kg s~ 2 (1 Jy = 10~ 26 W Hz -1 m~ 2 '). 

The luminosity distance £>l is defined in such a way that the conservation of energy applied to l(v) and s(v) takes 
the standard form of the continuity equation, 

J l{v)&v = A-kDI j s{u)Au. (Al) 

By definition a source is at redshift z, if electromagnetic radiation emitted by this source at a rest- frame frequency 
v c is observed at a frequency v = is e (1 + z)~ l \ or, in terms of wave-lengths, A Q = A (1 + z). A frequency interval dv 
around v c will be compressed to dv(l + when observed at v Q \ therefore, 



l{v e ) =47r J D 2 sK)(l + z)- i . (A2) 

If the source presents an emission line centered at a rest-frame frequency v c , one often considers the integrated 
luminosity and flux from the whole line. However, at least three definitions of these integrated quantities are commonly 
used. The most physically meaningful choices arc the frequency-integrated quantities, 

-\-Ais c (-^olAfo 

L = I l{v)dv, S= s(v)&v, (A3) 



v -Av 



where A^ e and /av q represent the half-widths of the line in rest-frame frequency and observer-frame frequency, re- 
spectively. The precise definition of these half-widths (i.e. the definition of where the line ends) depends on the 
observer's choice. L represents the actual power of the emission line and is measured in units proportional to 1W 
(1 L© = 3.839 • 10 26 W). S represents the power per unit area received by the observer, measured in units proportional 
to 1 Wm~ 2 = lkgs~ 3 . 

An alternative definition to the frequency-integrated quantities are the velocity-integrated analogues, often preferred 
by observers, 

n+AV f+AV 

L v = I l{v)dV, S v = / s(i/)dV, (A4) 



Ay 



Ay 



where V is the rest-frame velocity, projected on the line-of-sight, of the emitting material relative to the center of the 
observed galaxy, and AU is the maximal velocity (rotation+dispersion) of the emitting material. L v is measured in 
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units proportional to 1 kgm 3 s~ 3 (1 Jy kms^ 1 Mpc 2 = 9.521 • 10 21 kgm 3 s -3 ), and S v is measured in units proportional 
to lkgms -3 (Uykms -1 = 10~ 23 kgms -3 ). 

Confusion sometimes arises in the definition of the velocity V in Eqs. (|A4[) , since several definitions of velocity are 
commonly used in the context of emission and absorption lines (see Fig. [8]): (i) the standard recession velocity V op t{v) — 
c{v c —v)jv, traditionally used by optical astronomers; (ii) the variation V ra dio(z / ) = c{v c — v)/v c , sometimes employed 
by radio-astronomers; (iii) the "intrinsic rest-frame velocity" W, representing the rest-frame velocity, projected on the 
line-of-sight, of the emitting material relative to the center of the observed galaxy. For the investigation of emission 
(or absorption) lines at high redshift, it is critical to specify, which definition of the velocity V is used in the definition 
of velocity-integrated quantities like T v and S . The most natural choice, which we adopted in Eqs. (j A4|l . is V — W. 
This is the only choice, which makes T v an intrinsic property, that does not depend on the observer's distance. 

In the rest-frame of the observed galaxy, the center of the emission line is at the frequency v e and V is computed as 
V = c(v — v c )/v c . Hence, 

dV 

— — = A c if v is in the rest-frame. (A5) 

In the observer's frame, the center of the emission line is at the frequency v Q and and V is computed as V — c(v — v a )/v . 
Hence, 

dV 

— — = A Q if v is in the observer-frame. (A6) 
dv 
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Fig. 8. — Different velocity measures used in relation to emission lines from galaxies. 



Sometimes, line luminosities are defined with respect to the surface brightness temperature Tb, which is defined as 
the (frequency-dependent) temperature of a black-body with the physical size of the observed source and providing 
an identical flux density s(v a ). In radio astronomy the energy distribution of a black body is commonly approximated 
by the Rayleigh-Jeans law, i.e. the power radiated per unit of surface area, frequency, and solid angle equals u(y) — 
2v 2 kT^c~ 2 , where fcb is the Boltzmann constant. For isotropic line emission at rest-frame frequency v e this implies 
l{v c ) — AirD 2 pQ,u(v c ) — %-Kv 2 k\ 3 T'QC~ 2 D\tt, whe re D a — D^(l + z)~ 2 is the angular diameter distance and 17 is the 
solid angle subtended by the source. Using Eq. (|A2|) . we then obtain 



rp /„ x C 2 *{Vo)v c 2 (1+^) 3 ( , 7 , 

Tb( " o) = 2k h n • (A7) 

Tb is an intrinsic quantity, which does not change with redshift z, as can be seen from s(v )/Q cx (1 + z)~ 3 . Often 
the brightness temperature intensity / of an emission line source is defined as the velocity-integrated brightness 
temperature, 

r+AV 

1= / T B (^)dV, (A8) 

J-AV 

giving units proportional to lKms -1 (lKkms -1 = 10 3 Kms _1 ). Alternatively, observers sometimes define the 
intensity / as the velocity-integral of the "beam-diluted" brightness temperature T m b, which is smaller than Tb if the 
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sour ce does not cover the whole beam. We also note that some authors use the symbol I for fluxes, which we label S 
(e.g. I Weiss et al.ll200'ij). The brig htness temperature luminosity L T is defined as the product of the intensity and the 
source area (|Solomon et~aT]|1997l i. 

r+AV 

L T = D 2 A Q / T B (v)dV. (A9) 



J-AV 

This implies that L T is measured in units proportional to lKm 3 s _1 (lKkms _1 pc 2 = 9.521 • 10 35 Km 3 s _1 ). 

Basic relations 

From the conservation law of Eq. (|Alj) . we directly find the flux-to- luminosity relations 

L = AnDl S. (A10) 

The velocity-integrated flux S v can be expressed in terms of the frequency-integrated flux S by using the Jacobian 
of Eq. JM]), 

5 V = A 5. (All) 

Similarly, the velocity-integrated luminosities L v can be expressed in terms of the frequency-integrated luminosity L 
(i.e. intrinsic power of the emission line) by using the Jacobian of Eq. (|A5[) . 



L V = X C L, (A12) 

Finally, from Eqs. (TATl IMl lATOl [ATT]) . 

L T = (Snk^XlL. (A13) 

Using the four basic relations of Eqs. (|A10l[A"l~3|), we ca n express any of the quantities L, S, L v , S* v , and L 1 as a 
function of any other. For example, Eqs. (|A10| lAlT| IA12j) imply that 

L Y = {1 + z)- 1 AttDIS y . (A14) 
Note that this relation differs from Eq. (|A10|) by a redshift-factor. Eqs. (|A101 |AT2| imply that 



or, using typical observer-units, 



L=—Dlu S v , (A15) 
c 



, = 1.040 • 10" d — ^ -. (A16) 

L Q VMpc/ GHz Jykms" 1 y J 



This is equivalent to Eq. (1) in lSolomon fc Vanden Boutl (|2005h . Similarly, Eqs. (TATOl lATTl |AT3|) imply that 



LT = ^-^Dl(l + z )-*S v , (A17) 



or, in observer units, 



= 3.255- 10 7 

x (l + z)~ 3 



( D ^ 



Kkms^^c 2 VGHz/ V M P C 

S v 



(A18) 



Jy km s 1 

which is identical to Eq. (3) in lSolomon fc Vanden Boutl (|2005h . 

B. BACKGROUND OF THE CO/H 2 CONVERSION 

To-date, most estimates of molecular gas masses in galaxies rely on radio and (sub-) millimeter emission lines of trace 
molecules, especially on emission lines associated with the decay of rotational excitations of the CO molecule. 

It is not obvious that the CO-lines trace H2, and this method has indeed a long history of controversy. From local 
observations in the MW, it has become obvious that molecular gas resides in loosely connected giant (~ 10 — 100 pc) 
"clouds" , which are generally composed of hundreds of dense " clumps" (~ lpc), hosting even denser "cores" ( ~ 0.1 pc), 
where new stars are born (see e.g. the Orion Molecular Cloud. fMaddalc na et al.lll986l ; lTatematsu et al.lll993f) . CO-line 
emission cannot be used as mass tracer of individual clumps and cores, since they are optica lly thick to rotational CO- 
emission lines as can be inferred from the intensity-ratios between different rotational levels (jBinnev fc M crrificldl l 19981 
Chapter 8). However, if averaged over entire clouds or galaxies (typically 10 4 — 10 7 clouds), CO behaves as if it were 
optic ally thin, in a sense that individual clumps and cores do not significantly overlap (in space and frequency) (|Walll 
2006) , and hence on these large scales CO-line luminosities are expected to become suitable tracers of the molecular 
mass. Compelling empirical support for this conclusion was provided by the very tight correlation between the virial 
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Obreschkow et al. 



masses, estimate d from sizes an d velo city dispersions, and the CO(l-0)-luminosities of 273 molecular clouds in the 
MW analyzed bv lSolomon et al.l (|1987l ). 
To convert CO-line luminosities into H2-masses, it is common to define the X-factors as 

X, = (Bl) 

where J is the upper rotational transition J—yJ—1, Nn 2 is to column number-density of H2-molecules (here, we exclude 
Helium), and Ij is the brightness temperature intensity [see definition i n Eq. (|A8I)1 of the CO( J— > J— l )-emission line. 
Alternatively, the CO/H 2 conversion factors are sometimes defined as (jSolomon fc Vanden Boutll2005l ) 

a j = -j^-, (B2) 

where Lj is the brightness temperature luminosit y [see de finition in Eq. (|A9j) ] of the CO(J — » J— l)-emission line. 
Note that the defin itions of A^h 2 an d Mh 2 in Eqs. (|B1[ IB2j) do not include a Helium fraction, but some authors (e.g. 
iDownes et al.l [19931) include a Helium fraction of ~ 36% in Nu 2 and Mh 2 , which makes their values of Xj and aj 
1.36-times larger. 

Since Mh 2 = Dj^ilN^m^, where is the mass of a H2-molecule, and Lj = D\£llj for all J > 1, we find that 
the two conversion factors are related by 

a } = X 3 toh 2 , (B3) 

or, in typical observer units, 

"J = 1-6 -IP- 20 *,! 

MotKkms-ipc 2 )- 1 (Kkms^cm 2 )- 1 ' 1 ' 

From Eqs. (|A"l2l |A13l lB2l [B3|> it follows that 

**-3l*** (B5) 

where Lj and vj respectively denote the velocity- integrated luminosity and the rest-frame frequency of the CO(J-^ 
J— l)-cmission line, vj can be calculated as uj = J vco where vco — 115 GHz is the rest-frame frequency of the 
CO(l-0)-line. Eq. (|B5I) can then be expressed in typical observer units as 

M H2 = 313 J- 2 X 3 L co(j-+r-i) . . 

M Q lO^tKkms- 1 ™ 2 )- 1 ' Jykms^Mpc 2 ' [ ' 

Other mass-luminosity and mass-flux relations for H 2 commonly found in the standard literature can be derived from 
Eq. (|B6p and the basic relations in Section IA.2I 

Both the theoretical and the empirical determination of this conversion have a lon g history in radio astronomy, 
and are still consider ed hig hly challenging problems at the present day (see overviews in lMalonev fc Blacklll988l IWalll 
20071 iDickman et all 119861). Different methods to measur e aj (or Xj) were summarized by IDownes et al.l (|1995 



Arimoto et alj ( 19961 ). and ISolomon k. Vanden Boutl (|2005l) . The latter suggest that a sensible average value for the 
MW is ai — 3.4M©(Kkms _1 pc 2 ) -1 , or a% — 4.8M Q (Kkms _1 pc 2 ) -1 , if helium is included in the definition of a,\. 

C. LINE EMISSION OF THERMALIZED CO WITH FINITE OPTICAL DEPTH 

The rotational states of a diatomic molecule, such as CO, can be represented in the basis {|J, m)}, where J > is 
the angular quantum number and m s { — J, J} is the magnetic quantum number. In the absence of external fields, 
the energy only depends on J via Ej = h p vco J (J + l)/2, where Uco = H5 GHz and is the rest-frame frequency of 
the transition J = 1— >0. In local thermal equilibrium (LTE), the occupation probabilities of these energy-levels are 
therefore given by 



n3 = WT) exp 



hp vco J{J + 1) 
2fc b T 



(CI) 



where gj — 2 J + 1 are the degeneracies lifted by the quantum number m, T e is the excitation temperature (which is 
identical to the kinetic gas temperature in LTE-conditions), and Z{T) is the canonical partition function, which 
ensures the normalization condition J^nj = 1. The partition function is approximated to < 1% for all T > 10 K by 

z(Tc)= 2fcblc 2 
h p v co 3 

The interaction between a state | J, m) and the electromagnetic field only permits transitions simultaneously changing 
J by ±1 and m by -1,0,1. To determine the electromagnetic emission emerging from the three transitions \J, m) — > 
\J — l,m'), where m! 6 {m,m ± 1}, we require a measure of the rates of spontaneous emission from \ J,m), induced 
emission from [J, m), and absorption b y \ J — l,m'). These rates are effectively described by the Einstein coefficients 
Aj £?j.j_i, and Bj-\j (defined in lBinnev fc M crrificldl lT998l . Chapter 8), which can be calculated directly from 
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the interaction Hamiltonian between the rotational states and the electromagnetic-field. From consider ations of a gas 
in LTE it follows that these coefficients are related via Aj,j-i oc J a Sj,j-i and gj-i£?,j-i,j = gjSj.j-i. iRiegerl (|1974j ) 
showed that scales with J as 

and hence 

£?: 7_i oc and -Bt_i i oc . (C4) 

J J 1 2J+1 J M 2J-1 V y 

Following iBinnev fc Merrifieldl (|1998l ). the "source function" Zj, which is proportional to the power radiated per unit 
frequency from the transition (i.e. the sum of the power from all the transitions \ J,m) — > \J — X,m')) in an 

optically thick medium, is then given by 

h oc oc j- , (C5 

Bj-i,j-nj Bj,^! exp _ 1 

and hence the frequency-integrated power in a medium with arbitrary optical depth tj is given by 

Lj oc Jlj [1 - exp(-rj)] 

oc [1 - expt-Tj)] • f , ( C6 ) 

where 

tj(T c ) oc J -1 Bj_i :J -nj ,Bj,j-i) 

V ^co J 2 \ ( h p v C Q J \ (C7) 



oc exp 



^ sinh^ 



From Eq. (I A13|) the brightness temperature luminosity is given by Lj oc Lj J 3 . In the particular case of an optically 
thick medium (rj — > oo) and high temperatures (khT e 3> h p vco J), Eq. (|C6|) then implies that Lj is independent of 
J, which is indeed one of the essential properties of brightness temperature luminosities. 



